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Abstract 



We present two models of competitive learning, which are respectively interfacial 
(i.e., where the interiors of domains are unaffected by the rules of the model) and 
cooperative (i.e., where the bulk as well as the interface of an individual domain 
is governed by the rules of the model) learning. This learning is outcome-related, 
so that spatially and temporally local environments influence the conversion of a 
given site between one of two different types. We focus here on the behaviour of 
the models at coexistence, which yields new critical behaviour and the existence of 
a phase involving a novel type of coarsening which is oscillatory in nature. In the 
discussion we speculate on, among other issues, the likely behaviour of the models 
away from coexistence. 
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1 Introduction 



The notion of an outcome is familiar in the social sciences, but less so in the physical 
sciences. It is used in the former to describe the notion of a result following a course 
of action, especially in the context of game-theoretic applications [Q. In this paper we 
use this concept to devise models which, while motivated by ideas in the social sciences, 
are of interest to physicists both because they lead to new phenomena, such as that 
of oscillatory coarsening, as well as because, via their connection with known physical 
models, they are able to provide indications of new critical behaviour. 

The motivation of our model comes from the phenomenon of competitive learning; 
imagine the existence of two species, distracted (D) and receptive (R). We postulate that 
the D species is slower to learn than the R species, but that, on the other hand, the 
R species is more quick to forget; that is, a token R site is more swayed by the short- 
term successes of its neighbours to convert to their species. This could be a model of 
conservative versus radical societies, where the former adapt relatively slowly to change, 
but are more stable once changed, whereas the opposite applies to the latter. 

We shall define and investigate successively an interfacial version of our model in 
Sec. 0, and a cooperative version in Sec. 0. The emphasis will be put on the phase 
diagram of the model at coexistence, i.e., when neither of the species is favoured over the 
other; we will there examine the types of order present in the model, and the universal 
characteristics of the transition points between these phases. In Sec. |, we will discuss 
our results. 



2 Interfacial model 

2.1 Definition and generalities 

We here introduce the first of our two models (which we term the interfacial model) and 
set up definitions and notations. Individuals sit at the sites (or nodes) of a regular lattice 
with coordination number z. We shall most often consider the ci- dimensional hypercubic 
lattice, for which z = 2d. We represent the efficiency of the individual at site i as a 
(time-dependent) Ising spin variable: 

_ f +1 if « is of type R at time t, ,^ 
^ I —1 if z is of type D at time t. 

The efficiencies are dynamical variables, which are updated according to the following 
dynamical rules. 

• Step 1: Majority rule for site types 

In this step, we update the efficiencies via a zero-temperature (ferromagnetic) ma- 
jority rule. In other words, where a site is surrounded by a majority of its own 
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type, it stays the same, while if it is surrounded by a majority of the other type, it 
changes. In the absence of a majority, the site flips type with probability 1/2. 

More formally, this rule consists of ahgning the efficiency rji with the local field 
acting upon it, according to 



The local field. 



+1 if hi{t) > 0, 

±1 w.p. 1/2 ifhi{t) = 0, (2.2) 

-1 iihi{t)<Q. 



hi{t)^YlVj{t): (2.3) 

is defined as the sum of the efficiencies of the z neighbouring sites j of site i, and 
Ti is the associated time step. 

• Step 2: Association of outcomes with sites 

In this step, we assign to each site i an outcome cTj, which also takes Ising-like values 
+1 and —1, corresponding respectively to success or failure. The (time-dependent) 
outcomes cTj are random but correlated with the efficiencies r}i, namely 

iir„{t) ^+1, then a,{t + T,) = [^\ ^t!^^, 

if77,(i) = -l, then<7,(i + T2) = {+} "^-P- (2.4) 
^ ^ ^ ' [—1 w.p. 1 — P-, ^ ' 

with T2 being the associated time step. The performance parameters p± characterise 
the strength of correlation between efficiency and outcome for each species. 

This step is meant to model the performance of a particular site, according to its 
type. Thus according to our rules, if p+ > p_, the sites of type R are more likely to 
succeed at learning some new process (thus adapting more quickly to new patterns) 
than sites of type D. 

• Step 3: Conversion of sites according to performances 

In this step, we model the 'fickleness' of the sites where, as mentioned above, we 
assume that the R sites are more vulnerable to the performances of their neighbours 
than are the D sites. 

We denote by N^{t) the number of R-type neighbours of a given site i, and by 
N^{t) = z — -/V^(t) the number of its D-type neighbours. We also denote by I^{t) 
[respectively, /^"(t)] the number of these which have a positive outcome: 



_ + Vjjt) _ z + hjjt) /+(^) ^ + l + 



(2.5) 
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Then we postulate that the efficiencies are updated according to: 



if,,M = +la„d^<i|,the„,,(* + r3) = {;; ^ J-^^ 

■f ,,W = -1 and i|<||, then ,,(« + r,) = {+; ^ -:^_(2,6) 

with T3 being the associated time step. The efficiencies are not updated, i.e., they 
are left unchanged, if none of the conditions given in eq. (|2.6| ) is satisfied. The 
convertibility parameters e± characterise the fickleness of each species. 

The above rule is meant to hold only if the individual at site i is able to compare 
the outcomes of both types of its neighbours, namely for N^{t) and N~(t) 7^ 0, 
i.e., N^{t) 7^ and z. This is why we call this model the 'interfacial' model, since 
dynamical evolution occurs only for interfacial sites i whose neighbours are thus 
not all of the same type. Put another way, the interiors of clusters of any type do 
not evolve according to its rules. 

Then, if > and > S-, we note that even though the R sites are globally 
more successful than the D sites, they are more susceptible to local performance- 
based conversions than are the latter. 



So far we have only defined three elementary steps of the dynamics of the interfacial 
model, whose time scales, ri, T2, are arbitrary. We shall consider the realistic regime 
where the time scale over which individuals change type, i.e., Ti or t^, is much larger than 
the characteristic time scale of their activity over which their outcomes are updated, i.e., 
T2. Throughout the following, we thus set for simplicity 

r2 ^0, n = T3 = 1, (2.7) 

and we consider an ordered sequential dynamics, obtained by a regular 'sweeping' of the 
system. 



In the regime (|2.7| ) a given realisation of the outcomes (Jj{t) will enter the rules (|2.6| ) 
only once. Hence steps 2 and 3 can be recast as effective dynamical rules involving the 
efficiencies rii{t) and the associated local fields hi{t) alone. Obviously, these rules only 
hold for the sites i which are interfacial at time t, i.e., when hi{t) is different from its 
extremal values ±z. They are of the form: 



+1 


w.p. 


W- 




-1 


w.p. 


1 - 




+1 


w.p. 






-1 


w.p. 


1 - 


-w.{K{t)). 



\ir]i{t) = +1, then 7]i{t + 1) = 

ifr/,(t) = -l, thenr7,(t + l) = ^ 3 :''2'y/,\,,, (2.8) 



The effective transition probabilities w±{h) can be evaluated explicitly, by enumerating 
the 2^ possible realisations of the outcomes cxj of the sites j neighbouring site i, and 
weighting them appropriately. Consider for definiteness the example of the square lattice. 
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We have 2; = 4, so that the extremal values of the local field are h = ±4, while its values 
at interfacial sites are and ±2. The corresponding transition probabilities read: 



w+{+2) = l-e+p^l-pl), 

w_{0) = £_p+(l-p_)(2-p+-2p_ + 3p_p+), 
W-(-2) = s^p+{l-pt). (2.9) 



In order to predict the behaviour of the interfacial model for generic values of the 
performance parameters p± and the convertibility parameters e±, it is advantageous to 
examine in most detail the symmetric situation 

P = P+=P-, e = e+ = e^, (2.10) 

and then to investigate the infiuence of the two biasing fields 

H=p+-p_, B = £_-e+. (2.11) 

Indeed in the symmetric case ( |2.10| ), neither of the types R or D is favoured. In other 
words, the rules of the symmetric model are invariant under a global fiip of all the spins. 
The possible dynamical phase transitions between various kinds of collective behaviour 
are therefore expected to take place in this symmetric situation, in analogy with the 
phase transition in the ferromagnetic Ising model, which takes place at zero magnetic 
field (symmetric situation) and low enough temperature. We will discuss at a qualitative 
level the implications of this phase diagram for the generic situation of non-zero biasing 
fields B or H, in Sec. ^. 

This symmetric situation models the case when the two species, despite being identical 
in all their properties, are nevertheless distinguishable as being of two distinct types. 
The parameters p and e then control the intensity of surface noise. We attach the label 
'surface' to them because, as mentioned above, all discernible effects of these parameters 
are restricted to the interfaces separating clusters of individuals belonging to either of 
the two species. 

In the rest of this section we shall explore the behaviour of the symmetric interfacial 
model in the regime ( |2.7] ), via an approximate analytical treatment (known as the pair 
approximation) and then by means of numerical simulations. The basic quantities to be 
considered hereafter are the magnetisation M and the energy E. These quantities are 
defined for a finite sample of N individuals, i.e., N sites (or nodes) and Nz/2 bonds (or 
links), as 

M=l^T.V^, E=^y:i^-n.V,)- (2.12) 

« (y) 

In the following we shall usually consider the mean values (M) and (E), where the 
brackets represent an average over the random initial configuration of efficiencies {f]i{0)} 
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and over the stochastic dynamical rules (|2.2|) , ( |2.4| ), and ( p.6|) - i.e., the whole 'thermal 
history' of the system. The magnetisation (M) is the mean efficiency of an individual, 
while the energy {E) is the proportion of 'disparate' bonds (ij), such that rji ^ rjj. In 
the case of the cooperative model in three dimensions (to be presented later), we shall 
be led to consider also the dimensionless specific heat C. In analogy with equilibrium 
systems, this quantity is defined as the variance of the energy E per bond, namely 



(2.13) 

It will be advantageous to put the present model in perspective with the kinetic Ising 
model investigated in refs. 0, This is a dynamical model for Ising spins T]i{t) on the 
square lattice, defined by a stochastic rule of the form (|2.8| ), namely 



The transition rates W{hi{t)) depend only on the local field hi{t) defined in eq. ( |2.3|) . 
They assume the most general form compatible with spin-flip symmetry, i.e., W{h) + 
W{—h) = 1. This form involves two parameters (0 < pi,P2 < 1): 

W{+4) = p2, 

Wi+2) = pi, 

W{0) = 1/2, 

Wi-2) = 1-pi, 

W{-4) = 1-P2. (2.15) 

On the ferromagnetic side (1/2 < pi,p2 < 1), this model contains several special cases, 
including the Ising model with Glauber djTiamics, the (noisy) voter model, and the 
majority vote model. Introduced by de Oliveira et al. in the context of a general 
investigation of nonequilibrium spin models, it was subsequently treated by Drouffe and 
Godreche |^ who interpreted the parameters pi and p2 as two effective temperatures, 
linked respectively to an interfacial and a bulk noise. As our symmetric model only 
contains interfacial noise, it should correspond to the above model along the p2 = 1 line, 
with 1/2 < pi < 1 being a measure of the strength of interfacial noise. On this line, 
the model is paramagnetic for a large enough surface noise (1/2 < pi < 3/4), while 
it exhibits a frozen, i.e., totally ordered, ferromagnetic phase for a weak enough noise 
(3/4 < pi < 1). The transition point (pi = 3/4, p2 = 1) corresponds to the voter 
model. This model has been investigated extensively, both by mathematicians |^ and 
by physicists 0, and is known to be critical in two dimensions. 

2.2 Analytical approach: pair approximation 

The pair approximation is a particular case of the cluster method, an analytical approach 
proposed long ago [0] as a systematic improvement over mean-field theory to provide 
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approximate solutions of statistical mechanical models (see ref. for a review). Contrary 
to the standard mean-field approximation, or site approximation, the pair approximation 
has the advantage of taking into account correlations between pairs of neighbouring sites. 
A dynamical variant of the pair approximation had been introduced by Dickman in 
the framework of a surface-reaction model. More recently, de Oliveira et al. have also 
used a dynamical pair approximation to investigate their two-parameter kinetic Ising 
model §. 

We propose the following fully dynamical pair approximation, which leads to closed- 
form coupled evolution equations for the mean magnetisation M{t) and energy E{t) in 
our symmetric interfacial model; in the preceding, as well as in what follows, the angular 
brackets implying averaging have been taken to be implicit. 

The basic object of the pair approximation is the configuration of the dynamical 
variables attached to a bond. In the present case of a spatially homogeneous system, 
with binary site variables rji, an (oriented) bond {ij) can assume four configurations. We 
introduce the corresponding probabilities 

P{++) = Probjr^j = -|-1 and rjj = +1} = x, 

P{ ) = Prob{?7j = —1 and rjj = —1} = y, 

1 — X — y 



P(H — ) = Prob{?7i = -|-1 and rjj = —1} 
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P(-+) = Froh{r]i = -1 and rij = +1} = 1 — -. (2.16) 

The pair approximation consists of reducing any observable, i.e., the expectation of 
any function of the efficiencies rji, to a function of x and y, by systematically neglect- 
ing dynamical correlations between the efficiencies of any two sites which are not first 
neighbours, as well as any higher-order correlations. 

For instance, the probability law of the efficiency of any given site reads: 

1 + X — y 



P(+) = Fioh{r]i = +1} = P{++) + P{- 



2 



P(-) = Probfa = -1} = P(-+) + P(— ) = ^ ^ + ^ (2.17) 

The mean magnetisation M and energy E hence read: 

1-E + M 1-P-M,, 
M = x-y, E = l-x-y, i.e., x = , y = . (2.18) 

Thus, knowing x and y is equivalent to knowing M and E. 

As a second example, consider the conditional probability P(cr|r), defined as being 
the probability that rji = a, given that rjj = r for one of the neighbours j of site i. We 
have: 



p(+i+) = ^i±±i = p(-i+) = ^t±l = '---y 

^ P(+) l+x-y' ^ ^ ' P(+) 1 + x-y' 

p(+i-)^g±z)^^---^ P(H-) ^(--) 



(2.19) 



P(-) 1-x + y' ' ' ' P(-) 1-x + y' 



7 



The standard mean-field approximation corresponds to setting P{+\a) = P{+) = 
(l + M)/2 and P(-|cr) = P(-) = (l-M)/2 for a = ±1. The variables x and y are then 
related by (x - y f - 2(x + + 1 = 0, i.e., E = {1 - M^)/2. 

Within the dynamical pair approximation, coupled evolution equations for x{t) and 
y{t), or equivalently, for M{t) and E{t), can be derived by enumerating all the configu- 
rations of efficiencies attached to a given bond and its neighbourhood, weighting each of 
them with the appropriate probability, and by the appropriate transition rate for each 
step of the dynamics, as given in eqs. ( p.2|) , (p.4|), and ( p.6|) . 

We shall give just one example of how our calculations go. Consider a ++ bond {ij), 
and apply step 1 of the dynamics to site i at time t, say. The rule ( [2?^ ) involves the local 
field hi(t), or equivalently N^{t) = {z + hi{t))/2, defined in eq. ( p.5|) . Consider the z 
neighbours of site i. We know that the neighbour j is such that ?7j(t) = +1. The efficien- 
cies rjk of the {z — 1) other neighbours k are treated as independent random variables, 
each of them being distributed according to the conditional probabilities P{rik\+). The 
random variable Nl^{t) is thus distributed according to the binomial law 



FToh{N+ = n} 




{n = l,...,z). (2.20) 



.yy-i 

The probabilities of occurrence of the various cases of the rule ( |2.2| ) are thus known. 
This procedure leads to coupled first-order differential equations, of the form 

dM (dM\ (dM\ dE 



dt \ dt L \ dt ' dt 




(2.21) 



where the contributions of the various steps of the dynamics appear additively. The 
combinatorial analysis involved in the expressions of these contributions depends strongly 
on the lattice under consideration. 



• Two-dimensional case 

In the two-dimensional case, on the square lattice, the various contributions to the dif- 
ferential equations ( |2.21| ) read: 



(dM\ 2E^M 



3(1 - M^) - AE^ 

(5p2 - 5p - 2)(3 + M^)E^ - 4j9(5j9 - 9)(1 - M^)E 



\dt (1 -M2)2 

/dM\ _ 4:ep{l-p)E^M 
\df). ~ (1-M2)3 
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'dE\ _ 4ep{l-p)E/ 



- 2{2p - 1)(1 + 3M^)E^ + 2p{p + 3)(1 - M^)E^ 

- 3(p2 + p + 1)(1 - M^E + (p2 + p + 1)(1 - M^f^ . 



(2.22) 

Although these equations cannot be integrated exphcitly, the dynamical phase dia- 
gram within the pair approximation can be extracted from them as follows. The steady- 
state values of the magnetisation and of the energy are determined by equating the right- 
hand side of eqs. ( |2.22| ) to zero. The criterion for the presence of collective behaviour, 
i.e., long-range order, is the instability of the fixed-point solution such that M = 0. 

A linear stability analysis thus allows us to predict the existence of a transition line 
in the p-e plane, given by 

5 

''^^^ ^ Ap{i-p){i + 7p + 3p^y ^^-^^^ 

This line has two endpoints, pd = 0.47598 and Pc2 = 0.83060, corresponding to tran- 
sitions as a function of p, if we set e = 1. The intermediate phase {pi < p < p2 and 
ec{p) < e < 1) is paramagnetic, i.e., disordered, while the rest of the parameter space is 
frozen, i.e., totally ferromagnetically ordered. 

• Three-dimensional case 

In the three-dimensional case, on the cubic lattice, the coupled differential equations for 
M{t) and E{t) are much more lengthy than eqs. ( |2.22| ). We prefer not to write them 
down in full. Within this approximation, we predict the same phase diagram as for the 
square lattice, with 

292 

^"^^^ ~ 3p{l - p){32 + 1217p - 2243^2 + 2437p3 - 245^^) ' ^ ' ^ 

This line again has two endpoints, pd = 0.53047 and Pc2 = 0.86908, if we set e = 1 as 
before. We contrast the above with the predictions of standard mean-field theory which 
does not predict a phase transition for the symmetric interfacial model, in either two or 
three dimensions. 

The qualitative phase diagram of the interfacial model is shown in Fig. 1, together 
with that of the cooperative model investigated later on. Quantitative data concerning 
the location and the nature of the corresponding phase transitions are provided in Table 1. 



2.3 Numerical results 

We have performed a numerical investigation of the interfacial model throughout the 
p-e plane, both on the square and on the cubic lattice. The quantitative measurements 
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have been performed for e = 1, where the outcomes play a maximal role in the ordering 
behaviour of the model via step 3. The limiting situations p = and p = 1 both 
represent deterministic outcomes, where both species are restricted respectively to being 
total failures and total successes; in this case steps 2 and 3 are rendered essentially 
irrelevant, and the majority rule in step 1 lets the ordering proceed simply according to 
the types of the sites. 

We therefore discuss in the following the interfacial model on the square and on the 
cubic lattice, for £ = 1, as a function of p. In both cases we see clear evidence of a 
disordered, paramagnetic phase for pd < p < Pc2, between two frozen phases. This 
global picture is in qualitative agreement with the predictions of the pair approximation 
described in Sec. 



• Two-dimensional case 



We have already underlined the analogy between our interfacial model on the square 
lattice and the two-parameter model investigated in refs. 0, along the pi = 1 line, 
corresponding to interfacial noise alone. The strength of noise is measured by p2 in 
that model, while it is measured by p and e in our model. This suggests that the 
phase transitions at Pd and Pc2, which separate the disordered phase from both frozen 
phases, are critical points which belong to the universality class of the voter model: the 
intermediate disordered phase {pd < p < Pc2) corresponds to p2 < 3/4, while the two 
frozen phases (0 < p < pd and Pc2 < P < 1) correspond to p2 > 3/4. 

Fig. 2 shows snapshots of the dynamics of the model at times t = 8, t = 64, and 
t = 512, with a random initial condition. In anticipation of our numerical result Pc2 ~ 0.7 
(see below), we have chosen a value of the probability p = 0.7 so that our snapshots 
represent the time evolution of our system at (or very near) the critical point. Just 
as in ref. |^, the plots show a portion (of size 256^) of a square sample (of size 512^) 
with periodic boundary conditions. These plots bear a strong resemblance to those 
corresponding to the voter model (see Fig. 5 of ref. [Q), thus confirming our expectations. 

In order to provide quantitative confirmations of these observations, we have studied 
the decay of the mean energy E{t), starting with a random initial configuration. The 
two-dimensional voter model is known to exhibit unusual critical behaviour 0, this 
phenomenon is closely linked to the fact that ci = 2 is the marginal dimensionality for 
Brownian motion, which is known to be recurrent for d < 2 and transient for d > 2. In 
particular, the mean energy E{t) of the two-dimensional voter model falls off very slowly, 
as 

where the numerator 7r/2 is universal (i.e., it is independent of the microscopic details of 
the dynamics), whereas the time scale to is not. 

We have plotted in Fig. 3 the inverse energy 1/E{t) against \nt, for various values 
of p in the vicinity of the phase transition at Pc2 ~ 0.70, gleaned from our snapshots of 
the system. Each curve is an average over 200 independent samples of size 500^. The 
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graphs curve downwards in the disordered phase, where the energy converges to a non- 
zero equihbrium value E^qip), while they curve upwards in the frozen phases, where the 
energy falls off as E(t) ~ the latter is because of an underlying diffusive type of 

coarsening behaviour [^j. For p = 0.70, the data seem to become asymptotically parallel 
to the dashed line with slope 2/71 [see eq. ( p.25| )], shown on the plot. This observation 
confirms that our interfacial model belongs to the universality class of the voter model in 
a strong sense, i.e., including the prefactor of the law ( |2.25| ). Moreover, we thus obtain 
a rather accurate estimate for the transition point: Pc2 = 0.70 ± 0.01. A similar analysis 
yields Pd = 0.56 ± 0.01. The numerical values of both transition points are listed in 
Table 1. 

• Three-dimensional case 

In three dimensions, on the cubic lattice, we again find evidence of an intermediate 
disordered phase for pd < p < Pc2- Contrary to the two-dimensional case, we do not 
have any a priori knowledge of the phase transitions at pd and Pc2- In fact we do not 
know of any three-dimensional critical phenomenon driven by interfacial noise alone. It 
is nevertheless to be expected that the critical behaviour at these phase transitions cor- 
responds to a generic fixed-point behaviour, with finite, non-trivial values of the various 
critical exponents. This is in sharp contrast to the two-dimensional case, which belongs 
to the very special universality class of the voter model, with its well-known logarithmic 
behaviour ( p.25| ). In particular, we expect a non-trivial power-law decay of the mean 
energy at the transition points, of the form 

E{t) ~ r^. (2.26) 

This expectation is corroborated by our numerical results. Fig. 4 shows a log-log 
plot of the energy E{t), for various values of p in the vicinity of Pd- Each curve is 
an average over 50 independent samples of size 100'^. The two extremal dashed lines 
show the expected behaviour in the frozen and disordered phases, while the intermediate 
one corresponds to the critical law ( |2.26| ), with an exponent Q ~ 0.1. This picture is 
confirmed by a more refined analysis whose results are plotted in Fig. 5. This graph 
shows a log- log plot of the effective exponent Qcs(t) (defined as being the negative of the 
slope of the least-square fit to the data shown in Fig. 4, over the range t/2 < t' < 2t). 
The different kinds of behaviour corresponding to the frozen and disordered phases and 
to the critical point, appear more clearly. The data for p = 0.40, 0.42, and possibly 0.44, 
curve upwards and eventually tend toward the value 1/2, characteristic of the frozen 
phase. Conversely, the data for p = 0.48 and 0.50 curve downwards and eventually tend 
toward the value 0, characteristic of the disordered phase. The data for p = 0.45 and 
0.46 seem to converge to a non-trivial critical exponent Q. We thus obtain the estimates 
Pd = 0.45 ± 0.01 and = 0.10 ± 0.04. A similar analysis leads to Pc2 = 0.865 ± 0.005, 
with an exponent Q compatible with the above value. Although our numerical data do 
not allow us to rule out logarithmic behaviour of the mean energy, they definitely point 
towards critical behaviour characterised by a power law of the form ( |2.26D , with a finite 
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exponent f2. The numerical values of the transition points and of the exponent f2 are 
listed in Table 1. 

3 Cooperative model 
3.1 Definition 

One of the most interesting features of the model that we have investigated in Sec. § is its 
totally interfacial behaviour. However, from the viewpoint of the sociological behaviour 
that we are trying to model, it lacks an important feature to do with learning from the 
failures of neighbours when all the sites concerned are of the same type. Thus, consider 
a site of type D surrounded by others of its own kind in the region of phase space where 
> it is then plausible that a majority of the surrounding sites could fail at any 
given iteration. However, according to the rules of the interfacial model defined in the 
preceding section, the effect of this would be ephemeral; in step 3 the central site would 
be converted to a site of type R, but at the very next time step, the majority rule of 
step 1 would ensure that it was converted back to a D site. Thus there would be no 
long-term learning of the central site from the failures of sites of its own kind, to mirror 
the sociological phenomenon of 'learning from one's own mistakes'. 

To incorporate such long-term learning, we introduce a modification to our earlier 
model. In essence this ensures that the learning from failure in step 3 has a cooperative 
aspect, that is the central site as well as its neighbours learn from their failures and 
convert collectively to the other species. This in turn ensures that the majority rule 
of step 1 does not interfere with the conversion based on learning in step 3, since the 
central site and its neighbours are now all of the same kind, i.e., R, in the case of the 
above example. This rule (hard cooperative rule) is embodied in step 3a below. 

A modification of the above would be to say that only those sites which had failed, 
as well as the central site, would convert to the other species: then if a majority of 
surrounding sites fail, those sites as well as the central sites would convert at the next 
iteration to the other type. This rule (soft cooperative rule) is embodied in step 3b below. 
It is indeed a 'softer' version of step 3a; an alternative way of framing such a rule would 
be to have a stochastic formulation where a noise would control the conversion of sites 
in step 3, in proportion to the original number of failures. 

In contrast with the interfacial model studied so far, these variants of our cooperative 
model incorporate both surface noise and bulk noise. In this respect, the cooperative 
model is a generic dynamical model: in two dimensions, and at a qualitative level, it can 
be mapped onto the phase diagram of the two-parameter kinetic Ising model described 
above. It will turn out, however, that our cooperative model possesses a much greater di- 
versity of behaviour than the two-parameter model, since in addition to the conventional 
ferromagnetic and paramagnetic phases, it manifests a novel phase where the coarsening 
is oscillatory. 

In the following we focus on step 3a. We speculate that the soft version of step 3b 
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would only lead to a crossover behaviour from the pure interfacial model of Sec. 2, to 
the hard cooperative model defined by steps 1 to 3a. The rules discussed above are 
formulated respectively by: 

• Step 3a: Cooperative conversion of sites according to performance (hard) 

We add to step 3 by considering the case when a site i is surrounded by neighbours 
which are all of the same kind, either all R, i.e., N^{t) = z, or all D, i.e., Nf{t) = z, 
with the notation of eq. (|2.5|). Then we postulate that: 



if N+{t) = z and /+(t) < z/2, then r]i{t + n) = r]j{t + rg) = -1, 

if Nr{t) = z and /"(t) < z/2, then r]i{t + n) = T]j{t + rg) = +1, (3.1) 

where j denotes all the z neighbours of the central site i. This step of hard coop- 
erative conversion always involves the maximal number z + 1 oi individuals in a 
cluster made up of any central site and its neighbours. 

• Step 3b: Cooperative conversion of sites according to performance (soft) 

In this case we postulate that the cooperative rearrangement described by the 
rule ( prrp only applies to the central site i, and to the z — lf{t) neighbours j of 
site i who failed at time t. Hence this step of soft cooperative conversion involves 
a variable number z + 1 — lf{t) of individuals. 

In the following, we analyse this cooperative model with the hard rule (step 3a), along 
previous lines, and find that although it is more generic in some respects, in other ways 
it exhibits some very novel behaviour to do with the intermittent conversions of sites. 

We notice that, in the regime defined in eq. (|2l7|), steps 2 and 3a can be recast as 
effective dynamical rules involving the efficiencies rii{t) alone, in analogy with eq. (|2.8|) . 
These rules are valid only for sites i which are bulk sites at time t (i.e., when hi{t) = ±z), 
and are of the form 

if hi{t) = +z, then r]i{t + 1) = r]j{t + 1) = -1 w.p. 11, , . 

if hi{t) = -z, then r]i{t + 1) = rij{t + 1) = +1 w.p. H, ^ ' ^ 

where it is understood that the efficiencies are not updated in the complementary cases. 
The probability 11 can be evaluated by enumerating the 2^ possible realisations of the 
outcomes aj and weighting them appropriately. We thus obtain 

r(l-#(l + 3p) forrf = 2, 

I (1 -p)4(l + 4p+ 10p2) forrf = 3. ^ ' 

3.2 Analytical approach: pair approximation 



In order to analyse the cooperative model with the hard rule ( |3.1| ) within the pair ap- 
proximation, we have to determine the extra contribution of step 3a of the dynamics to 
the differential equations ( ^.211 ). The combinatorial analysis involved again depends on 



the lattice under consideration. 
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• Two-dimensional case 



In the case of the square lattice, the analysis is as follows. Consider for concreteness the 
first line of eq. ( p. 11 ). If the central spin is rji = +1, the move takes place with probability 
Ilp{+)P{+\+Y; it changes the total magnetisation by NAM = —10, while it leaves the 
energy unchanged. If the central spin is ?7j = — 1, the move takes place with probability 
np(— )P(+|— it changes the total magnetisation by NAM = —8 and the total energy 
by 2NAE = -4. 

We thus obtain after some algebra 



3 




+ 30(1 - M^E'^ - 5(1 - M2)3 
2(1 -p)3(3p + 1)^4(1 + 3M2) 



;i -M2)s 



(3.4) 



A linear stability analysis of the disordered solution M = of eqs. ( |2.21| ), including the 
contribution ( |3.4|) , leads to the following prediction. Within the pair approximation, the 
model has a unique transition at Pc = 0.86805, between a disordered phase at p < Pc 
and a ferromagnetically ordered one at p > pc- This prediction for the transition point 
is listed in Table 1. 



• Three-dimensional case 

The case of the cubic lattice is very similar to that of the square lattice. The pair 
approximation again predicts the existence of a unique phase transition, at pc = 0.87319. 



3.3 Numerical results 

For the cooperative model, updated in an ordered sequential way with rule 3a, either 
on the square lattice or on the cubic lattice, we see clear evidence of an intermediate 
disordered, paramagnetic phase for Po < P < Pc (Fig 1). The phase for pc < p < I 
is ferromagnetically ordered, while the phase for < p < Po exhibits a novel kind of 
dynamical order that we call 'oscillatory'. These observations are illustrated in the two- 
dimensional case in Figs. 6 and 7, which show snapshots of the model for p = 0.1 
(oscillatory phase) and p = 0.9 (ferromagnetic phase), respectively. 



3.3.1 Ferromagnetic phase transition 

The phase transition at p = pc might be expected to be in the universality class of the 
ferromagnetic Ising model, since it demarcates a disordered phase from a ferromagneti- 
cally ordered one. We have investigated this ferromagnetic phase transition in two and 
three dimensions, by means of numerical simulations. 
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• Two-dimensional case 



On the square lattice, the situation is very clearly seen from Fig. 8, which shows a plot 
of the mean magnetisation (M) against p. The samples used, of sizes 100^ and 200^, 
were large enough to prevent the system from flipping between one ordered phase and 
another during the simulations. We have also checked that the data shown in Fig. 8 are 
not affected by finite-size effects; indeed the data for the two system sizes mentioned 
above cannot be distinguished from one another on the plot. The dashed line shows a fit 
of the form (M) A^p — pc)^^^, with the well-known magnetisation exponent (3 = 1/8 of 
the two-dimensional Ising model. The quality of the plotted fit confirms our expectation: 
the ferromagnetic transition we observe does indeed belong to the universality class of 
the ferromagnetic Ising model. We also obtain an accurate evaluation of the transition 
point, Pc = 0.873 ± 0.002, i.e., = (7.4 ± 0.3)10^^^ which is listed in Table 1. 



• Three-dimensional case 

On the cubic lattice, the situation is less clear cut as a consequence of huge finite-size 
effects. An analysis of the mean magnetisation, analogous to that of Fig. 8, is totally 
inefficient, even for the purpose of locating the transition. We have instead considered 
the specific heat C, defined in eq. ( p.l3|) . Fig. 9 shows a plot of the specific heat C 
(multiplied by the transition probability 11), against p, for cubic samples of size L^, with 
various values of L. Each data point corresponds to an average over at least 100,000 
time steps. It is natural to consider the product 11 C, instead of C alone, because the 
specific heat is due to bulk fluctuations, whose driving force is the cooperative conversion 
embodied in step 3a; these take place, as mentioned before, at each time step with 
probability 11 (cf. eq. ( |3.3| )). We anticipate that the transition probability 11 will be very 
small throughout the ferromagnetic phase, just as in the two-dimensional case. 

The data for 11 C exhibit a peak, whose position, height and width depend on L. The 
critical exponents a and u can in principle be extracted from these data by means of the 



finite-size scaling law []TU 

C^L°/^f((p-pJL1/'^), (3.5) 

which implies that the height of the peak scales as L"/*^, while its relative position 
Pmax{L) — Pc with respect to the genuine transition point as well as its width, both 
scale as L"^/'^. Fig. 10 shows a plot of the peak position Pmax{L) against L'^/" . The fit 
shown as a dashed line incorporates a linear correction to scaling, of the form Praax{L) = 
Pc + L'^/^^A + B L^"^), and uses the recent estimates u = 0.630 and to = 0.80 |Tl|| cor- 
responding to the universality class of the three-dimensional Ising model. Reasonable 
agreement is found, although the amplitude for corrections to scaling is unexpectedly 
large. We thus obtain the estimate pc = 0.82 ± 0.01 for the critical ferromagnetic transi- 
tion, i.e.. He = (11 ± 2)10~^. We conclude therefore that the transition is most probably 
weakly Ising-like, by which we mean that corrections to scaling are very large, or else 
that the critical region is anomalously small. 
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3.3.2 Oscillatory phase transition 

For small values of the parameter p, i.e., values of the probability 11 close to unity, the 
cooperative system exhibits an oscillatory phase, already mentioned in the context of 
Fig. 7. 

This phase in unusual in several respects, and we were led to formulate its phe- 
nomenology as a result of various numerical investigations. First, since IT is close to 
unity, the net outcome of every time step is approximately to change all the efficiencies 
into their 'opposites'. This feature of the microscopic dynamics will be ignored hereafter 
by either only looking at the system at even times, or by considering quantities such as 
the skew magnetisation 

MskewW = (3.6) 

Second, and more importantly, we see that after a transient period, the system is filled 
with clusters which keep evolving forever. This perpetual motion is illustrated in Fig. 11, 
which shows a plot of the evolution of the skew magnetisation against t, for a square 
sample of size 100^, withp = 0.1 (hence 11 = 0.990711). The skew magnetisation exhibits 
slow, irregular oscillations, with a rather well-defined amplitude (roughly independent of 
size) and period (scaling roughly as L). This collective oscillatory behaviour remains 
noisy in the p — > limit, because of the presence of 'floppy' interfacial sites, whose local 
field vanishes. In order to investigate the transition at p = po between this unusual 
oscillatory phase and the usual disordered (paramagnetic) phase, we have used the order 
parameter (M^) = {M^^^^), which is clearly non-zero in the oscillatory phase. 

• Two-dimensional case 

Fig. 12 shows a plot of the order parameter (M^) against p, for square samples of various 
sizes. Each data point corresponds to an average over at least 100,000 time steps. The 
curves clearly exhibit a common point, corresponding to the transition. In the oscillatory 
phase, the order parameter has a non-zero thermodynamical value (M^)^ a; 0.16, with 
little or no dependence on p. It is discontinuous at the transition point, equaling {M^)^^ ^ 
0.06 at the transition point itself, and vanishing in the disordered phase. In the vicinity 
of the transition point, the order parameter is found to obey a finite-size scaling law of 
the form 

(m2)^<1.((p-P,)L1/-°), (3.7) 

illustrated in Fig. 13. The best data collapse is obtained for p^ = 0.136 ± 0.003 and 
l/i/<, = 0.45 ±0.10. 

• Three-dimensional ceise 

On the cubic lattice, the situation is qualitatively similar. The order parameter is still 
discontinuous, with a thermodynamical value around (M^)^ ^ 0.16 in the oscillatory 
phase. The values of the transition point, Po = 0.42 ± 0.005, of the associated exponent 
characterising the size dependence, l/ug — 0.8 ± 0.2, and of the order parameter at the 
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transition point, (M^)^^. ^ 0.0017 (with a large error bar), have been read off from Fig. 14. 
Each data point again corresponds to an average over at least 100,000 time steps. 

4 Discussion 

Starting from simple and intuitive ideas about the nature of experiential learning, we have 
devised two models involving competitive dynamics. The ideas are based on everyday 
experience: individuals who are quick to learn are often those who are also quick to adapt 
to changing circumstances (R type), while those who are slow to learn (D type) tend to be 
more conservative about what they know, adapting relatively slowly to changes in their 
environment. Our models embody these ideas, and emphasise the fact that adaptivity 
may occasionally lead to an overly quick change of behaviour (known colloquially as 
'knee-jerk reactions'), and that in this sense slow and steady learners can occasionally 
get the better of their quick counterparts. 

Both models have the same essential features: first, a zero-temperature majority rule 
involving species, whose effect is to 'convince' the central site to convert to (or stay 
the same as) the types of its neighbours. Next, outcomes for some process representing 
success or failure are assigned to each of these sites in proportion to their type, and 
last, based on these outcomes, the central site could decide to switch to another type 
or stay the same. This last step is rather subtle, as it involves the comparison of two 
ratios, I^{t)/N^{t) and I^{t)/N^{t), so that for example even one individual who is 
successful can have a greater impact on his/her neighbour than say two individuals who 
are successful in a trio of the same type. It is also in this last step that our interfacial 
and cooperative models differ, in that the latter allows for long-term changes even within 
a sea of the same species, depending on local outcomes. 

This diversity of behavioural modelling is, however, only one aspect of our work. The 
resulting rich and novel collective behaviour is what, for the most part, our paper is 
concerned with. We have focused on the symmetric situation, where none of the species 
is favoured. In particular, the two-dimensional case of our interfacial model turned out to 
be a physically motivated analogue of the kinetic Ising model discussed in a more abstract 
fashion in the literature |]^, ^]. In three dimensions, the coarsening behaviour of the same 
interfacial model has been characterised by a novel dynamical critical exponent Q ~ 0.10. 
The latter exponent can be expected to describe the relaxation of the energy of generic 
critical three-dimensional kinetic Ising models, whose rules involve only surface noise. 
The exponent fl certainly deserves to receive more attention in the future, in either 
three or higher dimensions. In the cooperative model, our phase diagram involves, in 
addition to regions of more conventional behaviour, a phase whose coarsening behaviour 
we term 'oscillatory'. This is the first example to our knowledge of such a non-stationary 
coarsening phenomenon, where domain walls (whose interiors flash between one species 
and another) appear to 'breathe', contracting and expanding forever in a more or less 
regular way. 

We now come back to the effect of the dimensionality d on the phase diagram of both 
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our models. This question is far from being trivial, as some steps of the dynamics are 
rather subtle, so that their net effect is hard to predict a priori. We can, however, make 
some educated guesses about the behaviour of our models. Starting with the interfacial 
model, we see that its one-dimensional version is frozen for any value of the parameters 
p and e. Indeed, the dynamics of step 1 are equivalent to the zero-temperature Glauber 
dynamics for the ferromagnetic Ising chain, while the only effect of the other steps is to 
reinforce the intensity of surface noise, and the mobility of the kinks between the growing 
ordered domains. This prediction is corroborated by numerical simulations, which clearly 
show a fall-off of the mean energy of the form {E) ^ A{p)t~^^'^, again because of an un- 
derlying diffusive coarsening behaviour 0. The amplitude A{p) is inversely proportional 
to the mobility of the point defects discussed above. We indeed find that A{p) is a maxi- 
mum at A(0) ^ ^(1) ~ (Stt)"^/^ = 0.1995 (the exact result for the Glauber-Ising chain), 
while its minimum lies around A{0.5) ~ 0.1751. As the dimensionality increases, the 
intermediate disordered phase is expected to open up as soon as c? > 1, to get larger and 
larger (we have indeed Pc2 — Pci ~ 0.14 for d = 2 and 0.41 for d = 3), and to invade the 
whole phase diagram as ^ oo. This increase of the disordered phase with dimension is 
not quantitatively reproduced by the pair approximation. The latter approach is however 
far more powerful than the standard mean-field theory, which does not predict any sort 
of critical point in the interfacial model. The situation of the cooperative model, with 
the hard rule of step 3a, is opposite to that of the interfacial model, and rather similar to 
that of standard equilibrium statistical-mechanical models. The one-dimensional model 
is always disordered, as its dynamics is equivalent to that of the Ising chain at high 
temperature. As dimension increases, the sizes po of the oscillatory phase and 1 — pc of 
the ferromagnetically ordered phase are observed to increase slowly. 

We now mention the many avenues we have left unexplored, speculating on the likely 
behaviour of some of them, based on our understanding of the phase diagram of our 
models in the symmetric case. First, we relegate to future investigations a more detailed 
investigation of the interfaces in our interfacial model. This issue has geometrical aspects, 
such as the surface tension and its critical behaviour, roughening phenomena, the cluster 
size distribution, and so on, as well as purely dynamical aspects, such as persistence and 
ageing. Another subject suitable for further development is the comparison of various 
kinds of dynamics. For a given model, i.e., given (stochastic) dynamical rules to be 
applied at every site, one can use either a parallel or a sequential updating procedure. 
Sequential dynamics can be ordered (the lattice is swept in a given ordered fashion), 
disordered (lattice sites are chosen at random), or more complex, such as the sublattice 
dynamics, which consists in successively sweeping the even and the odd sublattices. In 
the case of equilibrium statistical mechanical models, these features affect neither the 
phase diagram nor critical properties. In the present situation, changing the microscopic 
dynamical rules (inside a reasonable class) could have a mild but non-zero net effect on 
the phase diagram and on critical properties. 

Although we have restricted the investigations presented in this paper to the symmet- 
ric case when the performance parameters are identical, (because our main interest in the 
problem is largely to do with the critical behaviour at coexistence) the model is really 
meaningful from a sociological point of view only when p^. > p_, with the possibility 
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also of varying the convertibility parameters e_ and e+. In other words, we are led to 
consider the generic situation, where the biasing fields H and B, defined in eq. (|2.11| ), 



are non-zero. The case p+ > p_ and i.e., B < < H, corresponds to the R 

species being more 'genetically' successful, but more fickle than the D species, while for 
p+ > p- and i.e., both H and B positive, the R species is both more 'genet- 

ically' successful and more retentive than the D species. We would expect that in the 
latter case, the R species would usually be globally more successful in the long term, but 
that in the former case the final state of an initially random population may depend on 
the ratio H/\B\. There would be rather interesting conclusions from the point of view of 
societal evolution to be drawn from these investigations. 

Using the language of statistical mechanics, each of the biasing fields H and B is 
equivalent to a uniform magnetic field, favouring the R phase if it is positive. We expect 
that there is a non-trivial coexistence line in the [H, B) plane, above which the R phase 
is favoured (the case of sociological interest), and below which the D phase is favoured 
(its counterpart by symmetry). The fate of the D phase, whenever it is unfavoured, is 
dictated by the 'equation of state' of the model, relating the mean magnetisation M to 
the applied fields H and B. In the regime where both fields are small, the phase diagram 
of the symmetric situation investigated in this work allows us to make the following 
qualitative predictions: 

• In a disordered or paramagnetic phase, we have M ^ XiH + X2B in the ther- 
modynamic limit, with both susceptibilities being finite and positive. The mean 
population of the unfavoured species, P = (1 — M)/2, thus diminishes gradually, 
in proportion to the intensity of the biasing fields. In this regime, the coexistence 
line is linear: H = {x2/Xi)\B\. 

• In an ordered or ferromagnetic phase, there is a non-zero spontaneous magnetisation 
Mq. This means that the fraction of individuals in the unfavoured species, which 
is 1/2 in the symmetric case, decreases discontinuously to (1 — Mq)/2 whenever an 
'infinitesimal' biasing field is switched on. For a large but finite sample of linear size 
L, with N L'^ individuals, fields of order H B 1/N are sufficient to bring 
on the decay of the unfavoured species. The latter phenomenon is again driven by 
diffusion. It thus lasts a time of order r ~ ~ N'^^'^. 

• In a totally ordered or frozen phase, the spontaneous magnetisation assumes its 
maximal allowed value Mq = 1. In other words, all the individuals eventually 
belong to the favoured species, in the presence of an 'infinitesimal' biasing field. 
The finite-size estimates of the previous point still hold true. 

• Less can be said a priori in the case where the symmetric model is in the oscillatory 
phase, since this novel kind of order, where the growth of clusters is modulated by 
intermittent conversions, is not one where the usual phenomenology of equilibrium 
phase transitions can be readily invoked. We nevertheless speculate that the be- 
haviour of mean populations (possibly averaged over some intermediate time scale), 
will be qualitatively similar to that observed in a ferromagnetically ordered phase. 
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To close, we have presented in the above two models of competitive learning that 
are respectively interfacial and cooperative in their behaviour. We have focused on their 
properties in the symmetric situation where neither of the species is favoured over the 
other one, with emphasis on the phase diagram and on critical behaviour. While many 
sociological consequences of our models remain to be explored for the case when the 
species involved are unequally biased, we beheve that the behaviour at coexistence has 
already revealed a great deal of the underlying physics. The prediction of nontrivial 
critical behaviour entirely driven by surface noise in our interfacial model, or of the new 
phase of oscillatory coarsening in the case of our cooperative model, are possibly the most 
outstanding examples of this. We believe that further investigations of the issues sketched 
in the discussion could also yield phenomena, hitherto unknown, that are both novel in 
themselves and completely natural to the evolution of such truly complex systems. 
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Table and caption 



model 


transition points 
(pair approximation) 


transition points 
(numerical) 


characteristics 
of transitions 


2D interfacial 
(square lattice) 


Pel = 0.47598 
Pc2 = 0.83060 


Pel = 0.56 ±0.01 1 
Pc2 = 0.70±0.01J 


voter-like 


3D interfacial 
(cubic lattice) 


Pel = 0.53047 
Pe2 = 0.86908 


Pel = 0.45 ±0.01 1 
Pc2 = 0.865 ± 0.005 J 


n = 0.10 ±0.04 


2D cooperative 
(square lattice) 


Pc = 0.86805 


Po = 0.136 ±0.003 
Pe = 0.873 ± 0.002 


l/i^o = 0.45 ±0.10 
Ising-like 


3D cooperative 
(cubic lattice) 


Pc = 0.87319 


Po = 0.42 ± 0.005 
Pe = 0.82 ±0.01 


1/uo = 0.8 ±0.2 
weakly Ising-like 



Table 1: Characterisation of phase transitions shown in Fig. 1. Column 2: location of 
transition points within the pair approximation (see Sees. |2.2| and |3.2| ). Columns 3 and 
4: transition points and characteristics (exponents, universality classes) of transitions, 
obtained by means of numerical simulations (see Sees. ^]3| and p. 31) . 
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Figure captions 



Figure 1: Qualitative phase diagram of interfacial and cooperative models as a function 
of p. 

Figure 2: Snapshots of the dynamics of the interfacial model on the square lattice, for 
p = 0.7 ~ Pc2, at times (top) t = 8, (bottom left) t = 64, and (bottom right) t = 512. 

Figure 3: Plot of the inverse energy 1/E{t) of the interfacial model on the square lattice, 
against Int, for various values of p, indicated on the curves, in the vicinity of Pc2 ~ 0.70. 
The dashed straight line has slope 2/7i [see eq. ( |2.25| )]. 

Figure 4: Log-log plot of the energy E{t) of the interfacial model on the cubic lattice, 
against time t, for various values of p, indicated on the curves, in the vicinity of pd ~ 0.45. 
The dashed lines, meant as guides to the eye, have slopes (disordered), —1/2 (ordered), 
and —0.1 (critical). 

Figure 5: Log- log plot of the effective exponent ^les{t) of the energy of the interfacial 
model on the cubic lattice, against time t, for the same values of p as in Fig. 4. The 
symbol to the right with an error bar corresponds to the estimate f2 = 0.10 ± 0.04 for 
the critical exponent. 

Figure 6: Snapshots of the dynamics of the cooperative model on the square lattice, for 
p = 0.1, at times (top) t = 8, (bottom left) t = 64, and (bottom right) t = 512. 

Figure 7: Snapshots of the dynamics of the cooperative model on the square lattice, for 
p = 0.9, at times (top) t = 8, (bottom left) t = 64, and (bottom right) t = 512. 

Figure 8: Plot of the mean magnetisation (M) of the cooperative model on the square 
lattice, against p, in the ferromagnetic phase. The fit shown as a dashed line (see text) 
yields Pc = 0.873 ± 0.002. 

Figure 9: Plot of the specific heat C of the cooperative model on the cubic lattice 
(multiplied by the transition probability 11), against p near the ferromagnetic transition, 
for several sample sizes. 

Figure 10: Plot of the position Pmax(-^) of the peak in the data of Fig. 9, against L'^/" . 
Dashed line: fit discussed in the text, incorporating a linear correction to scaling. 

Figure 11: Plot of the skew magnetisation Ms^f.w{t)-, of the cooperative model on the 
square lattice, against time t, for a single history of a sample of size 100^, with p = 0.1. 

Figure 12: Plot of the order parameter (M^) of the cooperative model on the square 
lattice against p, for several sample sizes. 

Figure 13: Scaling plot of the data shown in Fig. 12, illustrating the finite-size scaling 
law (|3.7|) at the transition between the oscillatory and disordered phases, with po = 0.136 
and 1/uo = 0.45. 

Figure 14: Scaling plot of the order parameter (M^) of the cooperative model on the 
cubic lattice, illustrating the finite-size scaling law ( |3.7| ) at the oscillatory transition, with 
Po = 0.422 and l/uo = 0.8. 
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